Can the post-Newtonian gravitational waveform of an inspiraling binary be improved 
by solving the energy balance equation numerically? 
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The detection of gravitational waves from inspiraling compact binaries using matched filtering 
depends crucially on the availability of accurate template waveforms. We determine whether the 
accuracy of the templates' phasing can be improved by solving the post-Newtonian energy balance 
equation numerically, rather than (as is normally done) analytically within the post-Newtonian 
perturbative expansion. By specializing to the limit of a smaU mass ratio, we find evidence that 
there is no gain in accuracy. 



I. INTRODUCTION AND SUMMARY 

Several kilometer-scale interferometric gravitational 
wave detectors are currently being built, among them 
the two American LIGO detectors, the French- Italian 
VIRGO detector, the German-British GEO 600 detector 
and the Japanese TAMA detector. Gravitational waves 
from inspiraling compact binaries are among the most 
promising candidates to be detected. In order to extract 
a gravitational wave signal from the noisy background the 
technique of matched filtering will be used. One of 
the drawbacks of matched filtering is that the theoretical 
templates used must be close to the actual gravitational 
wave signal in order to detect the signal and estimate its 
parameters. In the case of a nearly circular inspiral of two 
point masses the expected gravitational wave signal has 
the form of a chirp, i.e., a roughly sinusoidal signal with 
gradually increasing amplitude and frequency. If such a 
signal is to be detected by matched filtering, high accu- 
racy is needed in the templates; in particular, the phase 
of the template must closely match the phase of the ac- 
tual signal. Inspiral templates have been calculated to 
date up to post-2. 5-Newtonian order [||. 

Several authors have investigated the question of to 
what post-Newtonian order does one need to push tem- 
plate computations in order to have an acceptably small 
template-inaccuracy reduction in event detection rate 
[Q-Q . The result is that post-2-Newtonian templates may 
be sufficiently accurate to detect neutron star/neutron 
star binaries the loss in event rate in this case is 
- 12% for initial LIGO and ~ 20% for advanced LIGO 

However, there are several motivations for trying to 
obtain more accurate templates. First, there may well 
be a high event rate of neutron star/black hole or black 
hole/black hole inspirals. For these more massive systems 
the accuracy requirements are more stringent, since the 
frequency band (~ 50 — 200 Hz) where most of the signal- 
to-noise is accumulated is in a more relativistic regime for 
more massive systems. For example, for initial LIGO de- 
tectors and for a binary system of a AMq black hole and 



a SOM© black hole, using post-2-Newtonian search tem- 
plates would allow us to detect only ~ 35% of the signals 
that otherwise would be detectable j^. A second moti- 
vation is that one will need high accuracy templates in 
order to avoid appreciable systematic errors in measure- 
ments of the binary's parameters pd]-p^. 

A variety of methods of increasing template accuracy 
have been pursued in recent years. First, one can com- 
pute the templates up to ever-higher post-Newtonian 
orders; this is arduous but going beyond the current 
post-2. 5-Newtonian templates is feasible. Progress is be- 
ing made on computing post-3-Newtonian templates; see 
Ref. [|lj] and references therein. Second, a celebrated 
result in this field was the discovery by Damour, Iyer, 
and Sathyaprakash [ pd] | that using Pade approximants 
can significantly increase the accuracy of template phas- 
ing. Third, Damour and Buonanno [ p5[ have suggested a 
particular ansatz for obtaining templates containing ad- 
ditional terms of all post-Newtonian orders, starting, say, 
from post-2-Newtonian templates. 

The purpose of this paper is to investigate yet an- 
other possible method of improving the accuracy of post- 
Newtonian templates. The basic idea is very simple. 
In computing post-2-Newtonian templates, for example, 
one should strictly speaking discard all terms of post- 
2. 5-Newtonian order (and higher) everywhere in the cal- 
culations. To do otherwise would be inconsistent with 
the perturbation expansion method. Yet, there could 
be pieces of the calculation for which retaining post-2. 5- 
Newtonian (and higher) order terms would lead to im- 
proved accuracy. For example, the dominant, m = 2 
piece of the waveform is usually written as 



h{t)^A{t) cos [0^^(t)] 



(1.1) 



where both the amplitude A{t) and the phase 
4P^ [t) have separate post-Newtonian expansions, A = 



Y.^e^A'^^^ and : 



(i) 



with e a formal expan- 



sion parameter. Now a perturbation theory purist would 
insist on inserting the expansion for (fP^ into Eq. (1.1) 
and on expanding the cosine using a Taylor expansion. 
However, it is well known that the resulting expression is 
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a much poorer representati on o f the true signal than the 
original un-expanded form ( |l.l| ). 

The question then arises: are there other stages in the 
construction of post-Newtonian templates where one dis- 
cards higher order terms, which, if retained, might lead to 
increased accuracy? A natural possibility is the stage in 
which one goes from the post-Newtonian formulae for the 
energy flux F{f) = —dE/dt{f) and orbital energy E{f) 
as functions of gravitational wave frequency /, to the 
formula for the phase (f)^^ (t) of the gravitational wave- 
form. Given analytical formulae for F{f) and E{f) up 
to some post-Newtonian order, one can either (I) solve 
analytically for within the post-Newtonian ap- 

proximation, discarding all higher order terms, or (II) 
one can numerically solve the energy balance equation to 
obtain This second procedure effectively gen- 

erates and retains terms at all post-Newtonian orders, 
so is strictly speaking inconsistent, but one might hope 
that it would lead to increased accuracy. We note that 
the papers investigating the accuracy of post- 

Newtonian templates have generally used method ( II) , 
whereas the popular data analysis package GRASP flfl] 
used in Ref. jl^l uses method (I). The GRASP manual 
p6[ speculates that method (II) might be more accurate. 

In this paper we present evidence, based on the limiting 
case of binaries with small mass ratios, that numerically 
solving the energy balance equation does not in fact in- 
crease the accuracy. We arrive at this conclusion after 
checking the accuracies of methods (I) and (II) in three 
different ways. We compare expansion coefficients of the 
Fourier transform of the waveform; we numerically find 
the relative error in the Fourier transform of the wave- 
form; and we compute overlaps of templates constructed 
from both methods with the exact waveform. The re- 
sult, that the numerical solution of the energy balance 
equation (method (II)) does not increase the accuracy, 
is disappointing, but constitutes useful information from 
the point of view of generating template banks for inspi- 
ral searches: there is no motivation in terms of increased 
event rate to solve numerically for the wave's phasing. 



II. METHOD OF CALCULATION 

In order to explain our calculation, we first summa- 
rize how the waveform's phasing can be computed from 
the energy flux function F{f) and orbital energy function 
E{f), where / is gravitational wave frequency pl],p|. Let 
mi , 777-2 be the masses of the two components of the bi- 
nary and 771 = 7711 + 77i2 be the total mass. Let (j){t) be 
the orbital phase of the binary, so that (t) — 2(j){t), 
where ip'~^^ is the phase of the dominant, m — 2 piece of 
the waveform. We define the dimensionless variable 



{■Kmf) 



1/3 



(2.1) 



[Here and throughout we use units with G = c = 1 .] The 
orbital phase (j){t) is derived from the relation 



and from the energy balance equation 
dE{v) 



dt 



-F{v). 



(2.2) 



(2.3) 



Equations ( |2.lD - (2.S) yield a parametric solution for 
(j){t) given by 



[V] ^ 0r / dvV 
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and 



t{v) = 



dv 



E'jv) 



(2.4) 



(2.5) 



where (j)^, tf. and Vi are constants. 

In the restricted post-Newtonian approximation, in 
which we neglect the m ^ 2 multipoles, the gravitational 
waveform has the form h{t) = A{t) cos[(j)'-^^ {t)] , where 
A{t) is a slowly varying amplitude. The Fourier trans- 
form h{f) of this waveform is h{f) = B{f)e^'^''^\ where 
B{f) is some frequency dependent prefactor, and where, 
in the stationary phase approximation, the phase ip{f) is 
given by 



7/;(/) = 27Tft{v) ~ 24>{v) 



(2.6) 



Using Eqs. (2^) and (2^) gives the frequency domain 
phase [||: 



iPif) = 2{t,/my - 24>c - 7r/4 

2 _ , 3 _3, E'{v) 



(2.7) 



We note that the corrections to the stationary phase 
approximation are very small, arising only at post-5- 
Newtonian order |l^, so it is sufficient for our purposes 
to use the ex pres sion (2.7). 

Equation (2.7) is the starting point for our analy- 
sis. We will investigate the accuracy with which vari- 
ous approximations reproduce the Fourier-domain phase 
?/'(/), which is the version of the phase function that 
is most relevant for matched filtering. The two possi- 
ble calculational methods we consider are (I) to insert 
post-Newtonian expressions for the functions E{v) and 
F{v) into Eq. ( |2.7| ), and discard all the higher order 
post-Newtonian terms generated, and (II) to insert post- 
Newtoni an e xpressions for the functions E{v) and Fiv) 
into Eq. (2/7) and solve exactly for the phase 7/'(/), retain- 
ing all the higher order post-Newtonian terms generated. 

To assess the accuracy of each of these two methods, 
we specialize to the limit 77717772/777^ — > for which the 
functions E{v) and F{v) are known P JTl| , pjt . We then 
check the accuracy of method (I) and (II) in three ways. 
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A. Checking the accuracy of methods (I) and (II) by 
comparing expansion coefficients of i/'(/) 

The first check is entirely analytical. We expand all the 
phase functions ip{f) as post-Newtonian power series in 
V up to some high order (e.g. post-5. 5-Newtonian), and 
compare the accuracy with which methods (I) and (11) 
reproduce the coefficients in this power series. While this 
comparison procedure is less accurate than comparing 
the phases produced by methods (1) and (II) to the exact 
numerical phase, it does allow us to check whether there 
is any indication that method (II) is more accurate than 
method (1). 

In more detail, our comparison procedure works as fol- 
lows. The post-Newtonian expansions for the functions 
E{v) and F{v) have the general form 



E{v) — ——rjmv^ 



E 



e2iV 



(2.8) 



Fiv) 



32 



2„,10 



-r] V 



i=2 



i=6 



(2.9) 



where rj = mim2l m? is the dimensionless mass ratio. 
The ellipses in Eq. ( |2.9|) represent possible terms propor- 
tional to (Inw)™ for TO > 2 which could arise at high 
post- New tonian ord ers. The coefficients e^, fi and gi in 
Eqs. (2.8) and (2.£) arc functions of the mass ratio rj. 
For general mass ratios, the coefficients and fi in are 
known up to 64 and /5J3], while for r] — al\ the co- 
efficients are known |^,^| and the fi and gi coefficients 
are known up to /n and g n [|l9|| . The known coefficients 
are tabulated in Appendix 



If we now insert the expansions (2_^) and ( |2.9| ) into the 
formula (2.7) for the phase ■(/'(/) we obtain 



128?7 



r., ^ 256r; o 128?/ . ^ 
dm 6 



(2.10) 



Here tx and K are constants which correspond to the 
initial time and initial phase, and the function P{v) has 
the expansion 

P{v) = l + ^[pj + qj\n{v) + rj[\n{v)f + . . .} i;^ 

(2.11) 



The coefficients Pj, qj and Vj in Eq. (2.11) are functions 
of the coefficients e^, fi and gi in Eqs. ( p.8[ ) and (2.E) for 
« < j: 

Pj =Pj{ei,---ejJi,...fj,gi,...gj), (2.12) 



1j ^qjiei,---ejji,...fj,gi,...gj), (2.13) 



Tj =rj{ei,...ejji,...fj,gi,...gj). (2.14) 
For example, the expressions for the first few pj 's are 
20(262 - /a) 



P2 = 



P3 



-4/3 



and 



Pi = 10(/| - 262/2 + 364 - fi 



(2.15) 



(2.16) 



(2.17) 



Suppose now that the functions E{f) and F{f) are 
known to post-A^-Newtonian order. Then the coefficients 
Ci, fi and gi are known for < i < 2A^. If we now fol- 
low the usual method (I) to generate the ph ase f unction 
'(/'(/), we obtain an expansion of the form (2.11) where 
the coefficients are given by 



^nP] 



and 



N' 3 



Pj(ei, 












■ e, , /i ,.../,, 51 , 



■ /li ■ • ■ fj^giy ■ 



■gj) J<27V, 
J > 2iV, 



(2.18) 



■gj) j<'2N, 
3 > 2A^, 



(2.19) 



■ gj) j<2N, 

J > 2N. 



(2.20) 



Here the superscript (I) means method (I) and the sub- 
script N refers to the post-A^-Newtonian approximation. 
On the other hand, if we use instead the method (II) to 
generate ip{f)y we obtain an expansion with expansion 
coefficients 



{II) 
nPj 



and 



Pj{ei,...ejji,...fj,gi,...g.j) j < 2N, 

Pj{ei, . . . 62JV, 0, . . .0, /i, . . . f2N, 

0,...0,5i,...g2JV,0,...0) j>2N, 

(2.21) 



qj(6i,...6j,/i,.../j-,gi,....gj) j < 2N, 

Qjiei, . . . e2N, 0, . . . 0, /i, . . . f2N, 

0,...0,gi,...52iv,0,...0) 3>2N, 



(2.22) 



rj(6i, . . . 6j, /i, . . . fj,gi, . . . gj) j < 2N, 
rj{ei, . . . 62JV, 0, . . . 0, /i, . . . /2JV, 

0,...0,5i,...52Ar,0,...0) j>2N. 



(2.23) 
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Thus, the two methods agree on pj, qj and Vj for j < 2iV, 
but for j > 2N method (I) gives expansion coefficients 
of zero while method (II) yields coefficients of the form 
Pj{ei, . . . e2N, 0, . . . 0, /i, . . . f2N, 0, . . . 0, 51, . . . g2N, 0, ... 0) 
which differ somewhat from the true values 
Pj (ei , . . . Cj , /i , . . . /j , gi , . . . (?j ) because of having the co- 
efficients Ci, fi and gi set to zero for 2N + 1 < i < j- We 
define 



_ (11) 



(2.24) 



and similarly for qj and rj, so that pk.k = Pk- 

As an example, suppose that the functions E(f) and 
F{f) were known only up to post-1.5-Newtonian order, 
so that only the coefficients 62, /2 and /s were known, 
but not 64 and f^. Up second post-Newtonian order the 
expansion (2.1f ) has the form 



P{v) = 1+P2V^ +P3V^ +P4V^ +0{V^), 



(2.25) 



where the coefficients P2, P3 and p^ are given in Eqs. 
(2.15) - ( ^.17 ) above. How accurately could we determine 
the coefficients p2 , Ps and p4 in this case? Obviously we 
could compute p2 and ps exactly since they do not depend 
on 64 and f^. However, the coefficient p^ does depend on 
64 and /4, and can be written as [cf. Eq. ( ^.17 ) above] 



P4 = P4,3 + Ap. 



'4,3- 



Here 



Pi,3 



{II) 
I.5P4 



10 (/I 



262/2 



(2.26) 



(2.27) 



is the piece of p4 that can be computed from the post- 1.5- 
Newtonian pieces of E{f) and F{f); it i s thu s nonlinear 
in the coefficients 62 and f2- The value (2.27) is the pre- 
diction of method (II), while the method (I) gives instead 
the value ^^5^4 = 0. The error term in Eq. ( |l26| ) IS given 
by 



Ap4 



10 (364 - U) 



(2.28) 



and is linear in the post-2-Newtonian coefficients 64 and 
/4. Using the values of 62, /2, 64 and listed in Ap- 
pendix^ we find that Ap4^3/p4 « —1.73 for ?7 = 0, which 
is rather large. Hence in this particular example we do 
not improve the accuracy in the coefficient p^ by using 
method (II) rather than method (I). 

In general, the question we want to address is whether 
the approximate coefficient ^^^Jpj = Pj,2N is typically 
significantly closer to the true coefficient pj than zero is 
to Pj, for j > 2N, i.e., whether 



\P].2N -Pj\ < , r , f ,^ 

— ^ (a few tens of percent) 

Pj 



(2.29) 



(and similarly for qj and rj). In Tables |[ || and III 
below we list the values of the true coefficients pj , qj and 
rj and also the approximate coefficients Pj,k, qj,k and Vj^k 



for various values of fc, c omp ut ed f rom t he va lues given in 
Appendix fusing Eqs. (U) and We list the 

analytic expressions for these approximate coefficients in 
Appendix Examination of Tables |[ || and III shows 
that there is no tendency for the inequality (2.29) to be 
satisfied. 

Therefore method (II) does not seem to lead to a gain 
in accuracy when compared to method (I) in the test 
mass case (?y ^ 0). 
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TABLE I. The coefficients pj^k = k/2Pj *s calculated according to method (11). These coefficients are what one ob- 
tains if the orbital energy E{f) and gravitational wave luminosity F{f) as functions of frequency / are known only up to 
post-fc/2-Newtonian order. Note that the values of pj,k differ significantly from their true values Pj = Pjj for k < j. 



method (II) 


k=2 


k=3 


k=4 


k=5 


k=6 


k=7 


k=8 


k=9 


k=10 


k=ll 


true values 


P3,k X 10"^ 





-0.0503 


















P3 X 10-3 


-0.0503 


P4,k X 10"^ 


0.0821 


0.0821 


0.0301 
















Pi X 10-3 


0.0301 


P5.k X 10"^ 





0.331 


0.331 


0.161 














PS X 10-3 


0.161 


P6,k X 10"^ 


-0.609 


-3.77 


-3.60 


-3.60 


-0.441 












P6 X 10-3 


-0.441 


P7,k X 10-3 





7.59 


7.52 


2.98 


2.98 


0.954 










P7 X 10-3 


0.954 


ps,k X 10-3 


-0.502 


-7.26 


-7.19 


-2.91 


0.828 


0.828 


0.995 








P8 X 10-3 


0.995 


P9,k X 10-3 





-37.8 


-40.2 


-28.8 


5.61 


11.6 


11.6 


4.43 






P9 X 10-3 


4.43 


PlO.fc X 10-3 


1.68 


43.3 


46.5 


15.2 


-1.76 


-12.0 


-11.5 


-11.5 


-8.77 




PlO X 10-3 


-8.77 


Pll,k X 10-3 





-76.8 


-82.5 


-28.4 


19.4 


26.1 


23.9 


14.4 


14.4 


12.3 


Pll X 10-3 


12.3 



TABLE II. The coefficients qj,k', see caption of Table ^. 



method (II) 


k=2 


k=3 


k=4 


k=5 


k=6 


k=7 


k=8 


k=9 


k=10 


k=ll 


true values 


g3,fc X 10-3 
























qs X 10-3 





q4,k X 10-3 

























g4 X 10-3 





gs.fe X 10-3 





0.992 


0.992 


0.482 














gs X 10-3 


0.482 


ge.fe X 10-3 














-0.326 












ge X 10-3 


-0.326 


g7,fe X 10-3 




























gr X 10-3 





g8,fc X 10-3 


1.51 


21.8 


21.6 


8.74 


-2.91 


-2.91 


-3.18 








g8 X 10-3 


-3.18 


g9,fe X 10-3 














-4.10 


-4.10 


-4.10 


-2.05 






gg X 10-3 


-2.05 


gio.fe X 10-3 














1.95 


1.95 


0.702 


0.702 


0.235 




gio X 10-3 


0.235 


gii.fe X 10-3 














-6.00 


-6.00 


-3.05 


-0.356 


-0.356 


-1.41 


gii X 10-3 


-1.41 



TABLE III. The coefficients rj^k', see caption of Table |. These coefficients vanish for j < 7. 



method (II) 


k=2 


k=3 


k=4 


k=5 


k=6 


k=7 


k=8 


k=9 


k=10 


k=ll 


true values 


rs.fc X 10-3 














1.29 


1.29 


0.584 








rs X 10-3 


0.584 


rg.k X 10-3 






























rg X 10-3 





rio,fc X 10-3 































rio X 10-3 





rii.fc X 10-3 
































ni X 10-3 
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B. Checking the accuracy of methods (I) and (II) 
numerically 

Next we perform a direct numerical check by compar- 
ing the phases produced by methods (I) and (II) to the 
exac t numerical phase. Note that the phase ip{f) in Eq. 
(2.7) is not directly observable, since it contains the un- 
known integration constants (j)c and tc, i.e. is de- 
termined only up to a linear function in /. Hence the 
relevant quantity to consider is 



2Tr^mE'{v) 
3w2 F{v) ' 



(2.30) 



In the test mass case E'{v) and F{v) are known exactly 
P, p]jis[| , and we can therefore find the exact 'ip"{f). 

Now suppose E'(v) and F{v) are only known up to 
post-fc/2-Newtonian order. If method (I) is used ip"{f) 
is found to be 



3w2 



E'{v) 
F(v) 



On the other hand method (II) yields 



[Fiv)]^ 



(2.31) 



(2.32) 



Here denotes the powerseries of the expression inside 
the brackets with terms kept up to order v'^. 

It is convenient to define the logarithmic relative errors 



In 



and 



i%r{f)-r{f) 



= In 



\E'(v)] 


E'iv) 


[f{v) \ 





F(v) 



(2.33) 



In 



iy>"(/)-^"(/) 



rif) 



= In 



E'iv) 
F(v) 



E'jv) 
F{v) 



(2.34) 



These errors are shown in Figures - It can be seen 
that there is no systematic tendency for method (II) to 
perform better than method (I). At post-2.5 and post-4- 
Newtonian order method (I) does better than method (II) 
for all V, while at post-3 and post-3. 5-Newtonian order 
method (II) is more accurate than method (I) . We would 
expect the same trend to hold for general values of the 
mass ratio rj. 




-0.8 -0.6 -0.4 -0.2 

(l/3)log(7tmf) 

FIG. 1. The errors in the phase tp{f) of the Fourier trans- 
formed waveform produced by methods (I) and (II) in the 
test mass limit. Plotted here are the logarithms of the rela- 
tive errors in the second derivative tf}"{f), for the case when 
the energy E{f) and gravitational wave luminosity F{f) are 
known only up to post-2. 5-Newtonian order. The horizontal 
axis is log(7rm/)/3 where m is the total mass of the system 
and / is gravitational wave frequency. The line denoted ISCO 
indicates the location of the innermost stable circular orbit. 
It can be seen that method (I) is more accurate for all fre- 
quencies /. 




-0.8 -0.6 -0.4 -0.2 

(l/3)log(7tmf) 

FIG. 2. The errors in the phase i/)(/) when E{f) and F{f) 
are known only up to post-3-Newtonian order; see caption of 
Fig. 1. In this case method (II) is overall more accurate. 
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-10 I ' • • ' • ' • ' 1 

-1.2 -1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 



(l/3)log(7tmf) 

FIG. 3. The errors in the phase V(/) when E{f) and F{f) 
are known only up to post-3. 5-Newtonian order; see caption 
of Fig. 1. In this case method (II) is more accurate for most 
frequencies /. 
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-10 I 1 • • 1 • 1 • 1 1 

-1.2 -1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 



(l/3)log(Jtmf) 

FIG. 4. The errors in the phase ?/'(/) when E{f) and F{f) 
are known only up to post-4-Newtonian order; see caption 
of Fig. 1. In this case method (I) is more accurate for all 
frequencies /. 




-1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 



(l/3)log(Jtmt) 



FIG. 5. The errors in the phase i>{f) when E{f) and F{f) 
are known only up to post-4. 5-Newtonian order; see caption 
of Fig. 1. In this case method (II) is more accurate for all 
frequencies /. 



C. Overlaps of templates constructed by methods 
(I) and (II) with the exact signal 

So far we have only considered how accurately the 
phase ip{f) is generated by methods (I) and (II). In this 
subsection we use the phases ^/2^(/) and ^/2^V'(/) to con- 
struct gravitational wave templates and then compute 
the templates' overlap with the exact waveform com- 
puted from the exact tp{f). 

We use the restricted post-Newtonian approximation 
and neglect the m ^ 2 multipoles. Thus the Fourier 
transform h{f) of the exact waveform is given by 

h{f) f-^/'^e'^'^^\ (2.35) 

Similarly we use methods (I) and (II) to construct the 
templates 

i%h{f) oc /-V6e^i%'A(/) ^2.36) 

and 

iy2^M/)oc/-V6e'iy>(/). (2.37) 

Next we compute Apostolatos' |2^] fitting factor (FF) 
to determine the templates' accuracy. The fitting factor 
is the ratio of the signal-to-noise ratio obtained with the 
imperfect template, to the signal-to-noise ratio that a 
perfect template would yield. The fitting factor can take 
values from zero to one, with unity indicating a perfect 
template. It is obtained from the ambiguity function 

A = ' ' (2.38) 

by maximizing over the template parameters, i.e. 

FF = max^^,t^A. (2.39) 

Notice that we hold the masses fixed in the maximization 
procedure: the templates and signal correspond to bina- 
ries of the same two masses. Here we have introduced 
the inner product 

JO 'Jri(/) 

where S'„(/) is the spectral density of the detector noise. 
The noise curve 5'„(/) used here is the Cutler-Flanagan 
fit [0 for the advanced LIGO. The largest contribution 
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to the overlaps comes from the frequency band between 
40 Hz and 200 Hz. 

We compute the fitting factors for several different 
choices of mi and m2 in order to get an indication of what 
might happen for general mass ratios, even though our 
results apply strictly only to the test mass limit f?-^ 0. 
The resulting fitting factors are listed in tables ^ and 
VI. Examination shows close agreement with the error 
plots of i^i^''^) and [/2V"(/) in Figs. |l| - |. At post- 
Newtonian orders where {j/2^"(/) has an error smaller 
than the error in ^y2^'0"(/) everywhere (e.g. post-2.5), 
method (I) always wins, and vice versa (e.g. post-4.5). 
On the other hand, at post-Newtonian orders where the 
error lines cross (e.g. post-3.5), the method with the 
smaller error in the w-region [v = (Trm/)^/'^] selected by 
the sensitive frequency band of the detector and the total 
mass m yields a larger fitting factor. 

Again we find that there is no systematic tendency 
for method (H) to be more accurate than method (I). 
Therefore method (II) does not lead to a gain in accuracy 
when compared to method (I). Our conclusion applies 
only to the limit — > 0, but we do not anticipate a 
different result for the general case. 

TABLE IV. The fitting factor (FF) at post-fc/2-Newtonian 
order for gravitational wave templates constructed by method 
(I) and (II). A fitting factor of unity indicates a perfect tem- 
plate. This table shows FF in the case of mi = m2 — 1.4Mq. 



k 


FF-method (I) 


FF-method (II) 


4 


0.495960 


0.710421 


5 


0.963230 


0.596752 


6 


0.985540 


0.981689 


7 


0.995678 


0.998342 


8 


0.998731 


0.997278 


9 


0.999264 


0.999851 


10 


0.999661 


0.999953 


11 


0.999911 


0.999964 



TABLE V. The fitting factor (FF) for method (I) and (II) 
in the case of mi — 1.4M0 and m2 — IOMq. 



k 


FF-method (I) 


FF-method (II) 


4 


0.329050 


0.871888 


5 


0.787947 


0.341547 


6 


0.826911 


0.872141 


7 


0.926867 


0.930399 


8 


0.967903 


0.919188 


9 


0.978698 


0.995905 


10 


0.980767 


0.990105 


11 


0.983868 


0.992731 



TABLE VI. The fitting factor (FF) for method (I) and (II) 
in the case of mi = m2 = IOMq. 



k 


FF-method (I) 


FF-nicthod (II) 


4 


0.483132 


0.908446 


5 


0.895806 


0.474364 


6 


0.896781 


0.976451 


7 


0.959834 


0.947991 


8 


0.968440 


0.948341 


9 


0.972900 


0.999764 


10 


0.975185 


0.994083 


11 


0.987979 


0.995887 
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APPENDIX A: COEFFICIENTS IN EXPANSIONS 
OF ENERGY AND ENERGY FLUX FUNCTIONS 

1. The coefficients a and /; up to 
post-2. 5-Newtonian order 



The coefficients in Eqs. ( p.8[ ) and ( |2.9j)_ up to 2.5 post- 
Newtonian order as given by BlanchetM are: 



62 



9 + 7? 

12 ' 



64 



27- 19?? + ryV3 



124 7 357? 
336 " T2"' 



/4 = - 



/3 = 47r, 



44711 92717? 65??2 



9072 



504 



18 



and 



„ /8 1 9 1 5357?\ 

^^-"l"67^ + ^j^- 



(Al) 

(A2) 

(A3) 
(A4) 

(A5) 
(A6) 
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2. The coefficients e^, fi and gi up to 
post-5. 5-Newtonian order 



The remaining coefBcients in Eqs. (2.8) and (2.9) have 
been given up to 5.5 post-Newtonian order in the test 
mass Umit in Ref. [H. These are 



66 



68 



610 



675 
""64 ' 

3969 

" 128 ' 

45927 
512 ' 



(A7) 
(A8) 
(A9) 



/6 = 



6643739519 1 7127 167r2 34241og(2) 



69854400 



105 



1628577 
504 ' 



105 



323105549467 2325977 136977^ 



/9 = 
flQ 



3178375200 4410 126 
39931 log(2) 47385 log(3) 
^ 294 1568 ' 

26597866751977 6848777 1369677 log(2) 



(AlO) 
(All) 

(A12) 
(A13) 



745113600 105 105 

2500861660823683 9166284677 42422377^ 



2831932303200 7858620 
83217611 log(2) 47385 log(3) 

11226600 ^ 196 ' 

839930975040177 177293777 

101708006400 ~^ 1176 
852128377log(2) 14215577 log(3) 



17640 

56 = 



784 



1712 



105 ' 

57 = 0, 

_ 232597 
~ 4410 ' 

684877 



and 



99 



gw 



911 



105 ' 

916628467 
7858620 ' 

17729377 
" 1176 



6804 
(A14) 



(A15) 

(A16) 
(A17) 
(A18) 

(A19) 

(A20) 

(A21) 



APPENDIX B: ANALYTIC EXPRESSIONS FOR 
THE COEFFICIENTS IN PHASE EXPANSION IN 
VARIOUS ORDERS OF APPROXIMATION 

1. Tlie coefficients Pj.n, Qj.n and rj,„ up to 
post-2. 5-Newtonian order for general mass ratios 

Here we Ust the analytic expressions for the coefficients 
Pj n, Qj n and rjn up to post-2. 5-Newtonian order for 
77^0. ■ 



_ 3715 5577 



P3,2 = 



P3,3 = -I677 



(Bl) 



(B2) 



(B3) 



5 (926521 -I- 1880368 7; -I- 905520 r;2) 
P4,2^P4,3=— ^^777^ — (B4) 



P-IA 



56448 



15293365 27145 7? 3085 77^ 
508032 ^ 504 ^ 72 



P5..2 = 



P5.,3 = P5A = 



20 (995 + 95277) 77 
189 



P5,5 



5 (7729 4- 252 7;) 77 
756 



92, n — q3,n — Q4ji — 



95,2 = 



95,3 = 95,4 



20 (995 -1- 952 77) 77 
63 



38645 77 ^ 



^2,n = r-s^n = r4 „ = Tg = 



(B5) 
(B6) 

(B7) 

(B8) 

(B9) 
(BIO) 

(Bll) 

(B12) 
(B13) 



9 



2. The coefficients Pj,n, qj.n and rj,n up to 
post-5. 5-Newtonian order in the test mass limit 

Here we list analytic expressions for the remaining co- 
efficients Pj^n, Qj,n and rj^n up to post-5. 5-Newtonian or- 
der in the test mass limit. 



P6,2 = -- 



5776858435 
9483264 



P6,3 



5776858435 
9483264 



320 TT^ 



37674179035 „„„ o 



P6,6 = 



1 0817850546611 6848 7 
4694215680 21~ 
640 TT^ _ 13696 log(2) 
3 21 



P7,2 = 



P7,3 



P7,i 



5680085 TT 
2352 



152000375 77 
63504 



P7,5 = P7,6 



241249475 77 
254016 



P7,7 



PS,2 = -- 



77096675 77 
254016 



7203742468445 
14338695168 



P8,3 = -- 



7203742468445 43160 tt^ 



14338695168 



63 



P&A 



P8,5 



499400855271485 43160 77^ 



1161434308608 



63 



499400855271485 47570 77^ 



1161434308608 



189 



(B14) 
(B15) 
(B16) 



(B17) 
(B18) 
(B19) 

(B20) 

(B21) 

(B22) 

(B23) 
(B24) 
(B25) 
(B26) 



Ps,6 = Psj 



3538122ir)!)41()7G17 1703440 7 



12775777394688 3969 
63110772 3406880 log(2) 



567 



3969 



(B27) 



P8,8 = 



2496799162103891233 36812 7 90490 77^ 

830425530654720 189 
1011020 log(2) 26325 log(3) 



567 



3969 



196 



P9,2 = 

-6756514105 77 , 

P9 3 = 640 77 

^^'^ 1185408 

-2308704875577 



P9A = 



P9,5 



3556224 
-971321608855 77 



- 64077^ 



640 77-^ 



P9,6 = 



341397504 
12113024196955177 27392777 640 77^ 



(B28) 
(B29) 
(B30) 

(B31) 

(B32) 



18776862720 
5478477 log(2) 
21 



21 



(B33) 



P9,7 = P9,8 = 



15706328988955177 27392 777 



18776862720 
640 77^ 5478477 log(2) 



21 



P9,9 



3 21 

90036665674763 77 13696 777 640 77^ 



(B34) 



18776862720 
27392 77 log(2) 

21 



21 



Pl0,3 



PlOA 



PW,5 



Pl0,6 



_ 1796613371630183 
~ 1070622572544 

1796613371630183 1240765 77^ 



+ 



1070622572544 294 

54094086068862461 11956093 77^ 

+ 



28906809458688 2646 

54094086068862461 458972531 77^ 
+ 



PW,7 



28906809458688 338688 

4027802547645341665 650561605 7 

317974904045568 333396 
312163997 77^ 650561605 log(2) 
435456 166698 

40278()2.')47()4.')341(i()5 650561605 7 

317974904045568 ^ 333396 
138083683 77^ 650561605 log(2) 



(B35) 

(B36) 
(B37) 
(B38) 
(B39) 



(B40) 



435456 



+ 



166698 



(B41) 



10 



Pl0,8 = Pl0,9 



23600127211067107843 116990189 7 



1878942614814720 



166698 



181984501 TT^ ^ 228376895 log(2) 



+ 



3048192 
15716025 log(3) 
21952 



333396 



(B42) 



Pio,io 



1412206995432957982751 6470582647 7 

+ 



126306697995878400 



27505170 



Pll,3 = 



578223115 TT^ 53992839431 log(2) 
^ 3048192 55010340 
5512455 log(3) 
21952 

Pll,2 = 

-40905234824185 TT 358720 tt^ 



Pll,4 



Pll,5 



7169347584 189 

-145661139175395577 _ 358720 tt^ 
193572384768 189 

-1211268636338065 71 112990 77^ 



Pll,6 = 



387144769536 189 
41090763354419749 7r 7577740 777 



Pn,7 



4258592464896 3969 
15130 7r3 15155480 7rlog(2) 
567 3969 

50239568645429749 n 7577740 7 ir 



+ 



+ 



4258592464896 3969 
15130 77^ 15155480 7r log(2) 



Pll,8 = 



567 3969 

3146788245124283189 TT 183628 777 

276808510218240 189 
94390 7r3 5572040 7r log(2) 

567 3969 
_ 26325 7rlog(3) 
49 



(B43) 
(B44) 

(B45) 
(B46) 
(B47) 

(B48) 

(B49) 



(B50) 



Pll,9 = P1I4O = 



1846304168796859019 TT 449308 777 



276808510218240 



3969 



94390 TT^ 1241720 TT log(2) 

567 ^ 3969 
26325 TT log(3) 
49 



(B51) 



Pii.ii = 



1795505143426433771 7r 3558011 777 
276808510218240 7938 
94390 7r^ 862549 7r log(2) 

567 1134 
26325 TT log(3) 
196 



(B52) 



96,2 = 96,3 = 96,4 = 96,5 = 



96,6 



6848 
21 



98,2 



97,n = 



7203742468445 
4779565056 



98,3 



98,4 



7203742468445 43160 77^ 



4779565056 



+ ■ 



21 



499400855271485 43160 tt^ 



387144769536 
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63 
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(B55) 



(B56) 



(B57) 
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